Thermal entanglement witness for materials with variable local spin lengths 
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We show that the thermal entanglement in a spin system using only magnetic susceptibility 
measurements is restricted to the insulator materials. We develop a generalization of the thermal 
entanglement witness that allows us to get information about the system entanglement with variable 
local spin lengths that can be used experimentally in conductor or insulator materials. As an 
application, we study thermal entanglement for the half-filled Hubbard model for linear, square and 
cubic clusters. We note that it is the itinerancy of electrons that favors the entanglement. Our 
results suggest a weak dependence between entanglement and external spin freedom degrees. 
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I. INTRODUCTION 

The manipulation of quantum systems in an entan- 
gled state that can be used as a quantum information 
channel is one of the main challenges of science today. 
Information theory, teleportation and cryptography are 
just some of the areas that may advance enormously 
through the amount of technological applications which 
can potentially make use of entanglement.— i 2 -^ It is also 
of great interest to explore the role played by entangle- 
ment systems in order to understand the basis of quan- 
tum mechanics.™^ However, even the quantification of 
entanglement remains an open question. Current re- 
searche is focusing on measures for precisely quantifying 
entanglement As an example, entanglement of indis- 
tinguishable particles calculated using different measures 
has shown that a same quantum state can have several 
different characterizations because of the lack of individ- 
ual identity of the entangled particles ^ 10 ' 11 . 

The use of uncertainty relations has provided an ef- 
ficient approach for obtaining one of the most precise 
experimental measures of entanglement^ 2 - A quantita- 
tive evaluation of the entangled states can be defined in 
terms of expectation values of a convenient witness oper- 
ator. This operator, called the entanglement witness, is 
defined as taking positive values for separable states and 
negative ones for entangled states. Thus, an appropri- 
ate uncertainty relation allows us to choose macroscopic 
properties which define an entanglement witness. 

A good level of interest has been focused on the spe- 
cial case of entanglement in macroscopic properties which 
has been particularly motivated by experiments that 
have shown the presence of entanglement in solid state 
systems . 13 ' 14 ' 15 ' 16 . Wiesniak et a/.— have recently ex- 
plored some aspects of the connection between entangle- 
ment and magnetic susceptibility for an arbitrary Hamil- 
tonian with spin length s. 

Experimental observations of thermal entanglement 
in spin systems using susceptibility measurements have 
been reported. Souza et al— have studied the com- 
pound Na2Cu5Si40i4. They found entanglement con- 
fined to the small clusters, with tripartite entanglement 



being stronger than bipartite entanglement. A similar re- 
sult was obtained by Vertesi and Bene 1 - 4 - in the Na2 V3O7 
system that formed a nanotubular structure of weakly 
coupled nine-site rings. Brukner et aZ.~ and Bose and 
Tribedi 16 showed entanglement in antiferromagnetic spin 
systems. 

Since these experimental susceptibility measurements 
are applicable to systems with spins localized in sites of 
the lattice, they are in accordance with the entangle- 
ment witness as defined by Wiesniak et a/.— However, 
this powerful tool is not adequate for systems with vari- 
able local spin lengths, which is an important feature for 
conductor materials. The present work addresses this is- 
sue. We will show that new aspects of the entanglement 
can appear when we consider systems with variable local 
spin lengths. The extension of the entanglement witness 
across variable local spins can be related to the itinerant 
electron models such as Hubbard and Falikov-Kimball. 

In this paper, we investigate the Hubbard models 
The purpose was to choose appropriated macroscopic 
variables in order to define an entanglement witness ad- 
equate for the Hubbard model, or any other model that 
can be applied to systems with variable local spin lengths. 

There are some works about entanglement associated 
to the Hubbard model^^Si 2 ^ 3 -^^ For example, 
the entanglement for the Hubbard dimer was investi- 
gated by Dowling et a/4 and Zanardi^. Gun et al. stud- 
ied the entanglement entropy on the extended Hubbard 
model and proposed that the entanglement can be used 
to identify quantum phase transitions. 2 ^ 1 Larsson and Jo- 
hannesson found exact expressions for the local entan- 
glement entropy on the one-dimensional Hubbard model 
at a quantum phase transition driven by a change in 
the magnetic field or chemical potential, related to the 
zero-temperature spin and charge susceptibilities. 21 Hu- 
dak modeled CeAl2 nanoparticles by the Hubbard model 
with negative chemical potential and, using entanglement 
entropy he studied the quantum phase transitions present 
in this system.— 

Some experimental results have indicated that the 
entanglement is restricted to small clusters within the 
materials ! 13 ' 14 Exploring this fact, we studied the critical 



2 



temperature below which there is thermal entanglement 
for finite chains and rings, using the standard direct diag- 
onalization method. 2 ^ 27 This approach is very well suited 
for small sized clusters since it produces exact results for 
thermodynamic quantities. Furthermore, it is also inter- 
esting to study the limit of large clusters. In this case, 
using the quantum Monte Carlo approac h 28 ' 29 ! 30 we ob- 
tained the temperature dependence of the entanglement 
witness for linear, square and simple cubic lattices as de- 
scribed by the Hubbard model. Summarizing, we will 
show how the cluster length, itinerancy of the electrons 
and system dimensions influence the thermal entangle- 
ment on the Hubbard model using direct diagonalization 
and quantum Monte Carlo methods. 

The organization of this paper is as follows. Entan- 
glement witness for constant and variable local spin is 
presented in Sec. II, the results in Sec. Ill, and the 
conclusions in Sec. IV. 



II. ENTANGLEMENT WITNESS 



The total magnetic susceptibility at null magnetic field 



X = Xx + Xy + Xz 



{^ B YNk B T 



(1) 



has been a useful variable to study the witness of thermal 
entanglement. Here, M is the total magnetization of N 
spins and (...) is the thermodynamic average. Consider- 
ing Si to be the length of the ith spin in the system, the 
entanglement condition for a thermal state of N spins of 
same length s (si = s for i — 1, N) is given byH 



X < 



knT 



(2) 



The above condition is deduced based on the method of 
entanglement detection using the uncertainty relations. 12 
In summary, an arbitrary thermal state of spin s has the 
follow condition 



Si 



8(8 + 1), 



< s 2 



(3) 



where Si is the spin vector of the individual site i. There- 
fore, if the thermal state is actually a product of N 
states of individual spins, the variance of magnetiza- 
tion would be the sum of variances of individual sites 



Si) > [s(s + 1) 



Nk B T X = Eti {S 
which is also valid for the general case of separable states 
due to the convexity of the mixture. 

However, we notice that the entanglement condition 
@ fails if the N individual spins have different lengths 
Si ^ s. Itinerant systems are an example of this phe- 
nomenon because the N individual sites can have differ- 
ent spin lengths due to the variety of ways in which they 
can be filled with particles. It can also occur in localized 



systems, since the sites can be filled in different ways. 
Thus, the Eq. © must be generalized as follows 



s, 



= N(L ) 



< 



(4) 



where s max is the largest spin length which the individual 
sites can take and Lj = -k Yli=i Si ■ Si + j is the spin 
spin correlation function. Therefore, we can rewrite the 
condition for entanglement ([2]) as 



X < 



{Lo)-8j 
k B T 



(5) 



Note that this is also valid for sites with same spin Si — 
s, because (Lq) = s(s + 1), s max — s and consequently 
the condition above is reduced to Eq. ©. 

Particularly, assume an TV-sites system in which the 
basis states are given by |nif, nii)(g>. . .<8>|njvr, 7i/vj.) where 
Tiia — or 1 (due to the Pauli exclusion principle) is 
the number of electrons with a-orientation of S z at the 
individual state i. Thus, Sj = for |0, 0) (vacuum state) 
or |1, 1} (singlet state of two electrons) and Sj = 1/2 for 
1,0) or |0, 1) (single electron states). Therefore, s max = 
1/2 and taking into account isotropy = L v - — L| (and 
consequently \ x — x y — X z )i the generalized condition 
of thermal entanglement can be expressed as 



£^X Z 



(LI) ~ 1/12 
knT 



< 0. 



(6) 



Note that if the individual state can only assume single 
electron states (|1, 0) or |0, 1)), Sj = s = 1/2 is fixed, 
(Lq) = 1/4 and the condition of entanglement above re- 
duces to © as hoped. 

The generalization of the entanglement witness for 
variable local spins introduces, besides the magnetic 
susceptibility, the Lq as an experimental measurement. 
Called local moment, the quantity Lq shows the degree of 
localization of electrons. This measurement is much less 
common and more difficult that the magnetic suscepti- 
bility. However, it can be obtained by neutron diffraction 
methods ^L 3 ^ 3 . 



III. RESULTS 

Using the witness (JSJ) , we investigated the thermal en- 
tanglement for the half filled itinerant electron systems 
described by the Hubbard model. The Hamiltonian is 



n = - l E ( c l«Cj* + hc) + U 



nifrni, 



(7) 



<ij>a 



where c\ a (ci a ) are the creation (annihilation) operators 
for electrons at site i, n,i a — c\ a Ci a , U is the on-site 
Coulomb (electron-electron) interaction and t is the near- 
est neighbor hopping integral representing the overlap of 
electron wave functions. 
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FIG. 1: Temperature dependence of the witness given by ex- 
pression ((6} for the four site linear chain of the half filled Hub- 
bard model using the grand canonical ensemble with U = 0, 
4 and 8. We adopt units ks = 1 and t = 1. 

We have obtained exact results for linear chains and 
rings with 2, 4 and 6 sites using the numerical method of 
direct diagonalization of small clusters over the canon- 
ical and the grand canonical ensembles . 26 ' 27 We have 
observed that the witness ^ for small odd numbers of 
sites provides no information about entanglement due to 
£ 3> for all T, since \ z diverges at null temperature.— 

It is illustrated in Fig. [T]that there is a critical tem- 
perature T c where £ (T c ) = and the system is entangled 
for T < T c , because £(T < T c ) < 0. Therefore, we can 
understand T c as the highest temperature below which 
the system is certainly entangled, since there is no cer- 
tainty about the entanglement when £ (T > T c ) > 0J£ 

Fig. [2] exhibits T c versus U for ID systems. A com- 
parison between results for different ensembles shows a 
good agreement at large Coulomb interaction, but not 
at small ones. We have found a rich dependence on U 
and N. With fixed N for small U/t, we see that the 
value of T c increases as the value of U/t is increased and 
tends to a maximum value. All curves have presented a 
value of interaction U which produces the maximum T c 
(global maximum of T c vs. U). For the strong Coulom- 
bian interaction U 3> t, we notice that all curves present 
a hyperbolic behavior T c oc U -1 . We will define the 
parameter rj(N,U) = A^kBT c j '{At 2 /U), where An is a 
function of N. r/(-/V, U) is convenient to compare our re- 
sults with the Heisenberg model ones, considering that in 
the asymptotic regime U 3> t there is an equivalence be- 
tween the half filled Hubbard and the Heisenberg models 
with exchange interaction J = At 2 /V \— 

For the grand canonical ensemble, as the size of an even 
sites system increases, the values of the maximum global 
l/ max and T i ™ ax also increase. Although the canonical 
ensemble has a similar increasing relation between {/ max 
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FIG. 2: Coulombian interaction dependence of the critical 
temperature for the finite one- dimensional half filled Hubbard 
model. We adopt units ks = 1 and t = 1. The solid and 
dashed lines are related to chains and rings, respectively. Each 
curve is labeled by its number of sites. 



and N, there is no monotonic behavior of T c max versus 
N. We performed a numerical extrapolation using the 
grand canonical ensemble for linear chains with 2, 4 and 
6 sites. Our extrapolation analysis predicts A j BT ( ! nax = 
0.712£ at f7 max = A.lt in the thermodynamic limit. We 
also obtained 77(00,00) = 1.568 ± 0.003 which is very 
close to the exact value ij(qo, 00) = ksT c / J = 1.6 for i-s 
Heisenberg models 

The direct diagonalization approach is very suitable 
for small sized clusters, but becomes inefficient when the 
system has its size increased. On the other hand, the 
quantum Monte Carlo (QMC) metho d 28 ' 29 is an efficient 
approach to study large systems. Using it we have stud- 
ied the entanglement witness for linear, square and cubic 
lattices. 

The QMC method treats the exponentials of the grand 
partition function with the Suzuki- Trotter decomposition 
scheme. Using a discrete Hubbard-Stratonovich transfor- 
mation it converts the electron-electron interaction into 
one of free electrons interacting with a time-dependent 
Ising field. With it we compute the relative weights of 
the Ising field configurations. The algorithm follows the 
lines of those for classical systems, except for the Boltz- 
mann weight that is expressed as a sum over Ising spins 
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FIG. 3: Coulombian interaction dependence of the critical 
temperature for the one, two and three dimensional half filled 
Hubbard model. We adopt units Ub = 1 and t = 1. The 
scatter is related to the quantum Monte Carlo method and 
the solid line is an extrapolation of the thermodynamic limit 
for a linear chain. 

of a product of determinants. For the Hubbard model at 
half filling, the product of determinants is always posi- 
tive. We have used the imaginary time discretization of 
the QMC At = 0.1253 

Fig. [3] shows the QMC results for a 64-site ring, and 
for 100-site square and cubic lattices. Our T c cannot 
be estimated accurately for large U/t because the QMC 
becomes unstable at low temperatures and with strong 
Coulombian interaction^ 

Furthermore, we include in Fig. [3] the thermodynamic 
limit extrapolation obtained through the linear chain re- 
sults from the direct diagonalization for small clusters. 
Note that the extrapolation is consistent with the simu- 
lation. Notice that for different lattices, the results for 
T c are similar revealing no new behavior. 

IV. CONCLUSION 

Wesniak et al<^ have suggested that magnetic suscep- 
tibility can be a macroscopic (thermodynamical) spin en- 



tanglement witness without complete knowledge of the 
specific model (Hamiltonian) of the solid. However, we 
observed here that its applicability is restricted to the in- 
sulator materials because local features of the spin length 
affect the deviation of the witness. We have developed 
a generalization that allows us to get information about 
system entanglement with variable local spin lengths such 
as found in itinerant electron systems. Moreover, our 
witness is also valid for fixed local spin lengths and con- 
sequently, it can be used experimentally in conductor or 
insulator materials. 

As an application, we studied thermal entanglement 
for the one, two and three dimensional half filled Hubbard 
model. We obtained the critical temperature T c below 
which the system is certainly entangled. We have shown 
that there is a Coulombian repulsion that presents a 
global T c maximum. This feature is relevant for quantum 
information science, since it reveals the optimal Coulom- 
bian repulsion referent to the highest temperature where 
the system is definitely entangled. In addition, the de- 
crease of T c for t <C U indicates that the itinerancy of 
electrons favors the entanglement. Furthermore, at the 
asymptotic regime U 3> t we show, through a numeri- 
cal extrapolation to the thermodynamic limit, that T c is 
in accordance with the exact result for the ^-s Heisen- 
berg model. A recent study has shown that higher spin 
length increases the T c rH Since higher spin length means 
higher internal degrees of freedom the above result shows 
a strong favoring of entanglement according to the in- 
creases in the internal degrees of freedom. In this work, 
from the results of T c for linear, square and cubic lattices, 
we notice that an increase in the external spin degrees of 
freedom produces similar results. These results suggest 
that the dependence between entanglement and internal 
spin degrees of freedom is far stronger than between en- 
tanglement and external spin degrees of freedom. 
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